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Context. We have studied the evolution of molecular gas during the early stages of protostellar collapse, when the freeze-out 
of 'heavy' species on to grains occurs. 

Aims. In addition to studying the freeze-out of 'heavy' species on to grains, we wished to compute the variation of the population 
densities of the different nuclear spin states of 'tracer' molecular ions, such as H2D + and D2H + , which are currently observed 
, only in their ortho and para forms, respectively. 
Methods. Chemical processes which determine the relative populations of the nuclear spin states of molecules and molecular 
ions were included explicitly. Nuclear spin-changing reactions have received much less attention in the literature than those 
leading to deuteration; but, in fact, the former processes are as significant as the latter and often involve the same reactants. A 
Q " 'free-fall' model of gravitational collapse was adopted. 

S—j | Results. The ortho:para ratios of some species, e.g. H2D + , vary considerably as the density increases. Because the dynamical 
C/5 timescale is relatively short, there can be large departures of the predictions of the free-fall model from the steady-state solution 
' at the same density. Values of the ortho:para H2 ratio much higher than in steady state, which would prevail in 'young' molecular 
clouds, are found to be inconsistent with high levels of deuteration of the gas. The internal energy of ortho-fb acts as a reservoir 
V of chemical energy which inhibits the deuteration of H3" and hence of other species, such as N2H + and NH3. 

' Conclusions. The principal conclusion is that the degree of deuteration of molecular ions and molecules is sensitive to the 
^ , ortho:para H2 ratio and hence to the chemical and thermal history of the precursor molecular cloud. 
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1. Introduction 

Studies of the early stages of protostellar collapse are 
presently at the leading edge of research in astronomy. 
These studies have been driven mainly by new observa- 
tions, which suggest that, during an initial, approximately 
isothermal phase of the collapse process, molecular species 
containing the heavy elements, notably carbon and oxy- 
gen, are removed from the gas phase by 'freeze-out' on 
to the surfaces of grains (Bacmann et al. 2002, Tafalla 
et al. 2004). Nitrogen-bearing species, specifically N2H" 1 " 
and NH3, subsist longer (Bergin & Langer 1997, Belloche 
& Andre 2004) but eventually freeze out also. There re- 
mains a gas which is 'fully depleted', in the sense that only 
hydrogen bearing species and helium remain. 

In conditions under which CO is heavily depleted from 
the gas phase, H + and H^" become the principal ions. The 
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abundance of H + is controlled by its recombination with 
electrons on the surfaces of negatively charged grains: two- 
body recombination in the gas phase is a radiative process 
and hence slow. Indeed, the relative abundance of H + and 
H3" during collapse depends on the surface area of neg- 
atively charged grains, per unit volume of gas, which is 
available for the recombination of H + ions. If the rate of 
coagulation of the grains is insufficiently rapid, compared 
with the rate of collapse of the protostellar core, the grain 
surface area remains high and the H + ions are more effec- 
tively neutralized; H^" becomes the dominant ion. 

Observations have shown, and calculations have con- 
firmed, that the abundances of the singly and multiply 
deuterated forms of Hjj" can attain exceptionally high val- 
ues under conditions of extreme heavy element depletion 
(Caselli et al. 2003, Ceccarelli et al. 2004, Vastel et al. 
2004; Walmsley et al. 2004, Ceccarelli & Dominik 2005). 
Furthermore, these high degrees of deuteration can prop- 
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agate to the residual nitrogen-bearing species, such as 
N 2 H + and NH 3 . The deuterated isotopes of NH 3 have 
been detected with fractional abundances which are many 
orders of magnitude higher than would be anticipated on 
the basis of the D:H elemental abundance ratio, which is 
of the order of 10 -5 . The related reaction networks have 
been considered in detail by several authors (Roberts et 
al. 2003, 2004; Roueff et al. 2005). Deuteration occurs in 
dcuteron-proton exchange reactions with the deuterated 
forms of H^. The analogous proton exchange reactions 
modify the relative population densities of the different 
possible nuclear spin states of molecules and molecular 
ions, such as NH 3 and H^ itself. These processes have 
received much less attention in the literature but are, in 
fact, as significant as deuteration. Observations of ortho 
and para forms could provide information on the physi- 
cal conditions in the core during collapse. Moreover, when 
only one of the possible nuclear spin states of a species is 
observed, which is usually the case, calculations of the rel- 
ative population densities provide a means of estimating 
the total number density of the species. 

The ortho:para ratio of species such as H^ and H 2 D + 
are dependent on the ortho:para H2 ratio and vary with 
the gas density in the course of the gravitational col- 
lapse. Furthermore, the degree of deuteration of Hj de- 
pends on the density of ortho-H2, through the reaction 
ortho-H 2 D + (ortho-H 2 , HD)!!^, which is the reverse of 
the process of deuteration of H^ (cf. Pineau des Forets 
et al. 1991, Gerlich et al. 2002). The analogous reac- 
tions involving para-H 2 and para-H 2 D + arc endother- 
mic and slow at low gas temperatures. However, para- 
H 2 D + is converted into ortho-H 2 D + in the reaction para- 
H 2 D+(ortho-H 2 , para-H 2 )ortho-H 2 D + , which is exother- 
mic. Thus, as the assumed value of the initial abundance 
of ortho-H 2 (equivalently, the initial ortho:para H 2 ratio) 
increases, both ortho- and para-H 2 D + are converted back 
to H^. Consequently, the initial degrees of deuteration, 
not only of Hjj~, but also of nitrogen-bearing species such 
as N 2 H + and NH 3 , fall to low values, as H 2 D+ is the main 
agent of their deuteration. It should be clear from these 
comments that the initial value of the ortho:para H 2 ratio 
is critical to the composition of the collapsing gas, partic- 
ularly to its level of deuteration. Indeed, the high levels of 
deuteration which are observed in some protostellar cores 
impose constraints on the initial ortho:para H 2 ratio. 

Molecular hydrogen forms on grains, at a rate which 
is still not well determined. It is believed (and assumed 
here) that the ortho:para ratio of the hydrogen molecules 
which are formed and released into the gas phase is the 
statistical value of 3:1. Subsequent proton-exchange reac- 
tions with H + and interconvert ortho- and para-H 2 . 
However, at low temperatures, the timescales associated 
with these reactions are large, owing to the low degree 
of ionization of the gas; this is particularly true of para 
to ortho conversion, which is endothermic. As we shall 
see below, it is certainly not clear that the ortho:para H 2 
ratio at the commencement of gravitational collapse of 
the molecular gas has attained its value in steady state; 



and it is almost certainly higher than its value in thermal 
(Boltzmann) equilibrium. At low temperatures, ortho-H 2 
is formed principally by H-atom recombination on the 
surfaces of grains, rather than through proton-exchange 
reactions with para-H 2 in the gas phase. Accordingly, we 
have treated the initial ortho:para H 2 ratio as a parameter, 
which we have varied between its upper limit of 3 and its 
lower limit, taken to be the ratio of the populations of the 
J = 1 and J = rotational states in thermal equilibrium 
at kinetic temperature T. A high fractional abundance of 
ortho-H 2 constitutes a major source of (internal) energy 
in molecular gas at low kinetic temperatures (T w 10 K): 
the lowest state of ortho-H 2 , J = 1, lies approximately 
170 K above the ground state of para-H 2 , J = 0; this en- 
ergy becomes available in proton-exchange reactions with 
ortho-H 2 . 

In Section 2, we outline the model of gravitational col- 
lapse which we have employed and consider the processes 
which influence the ortho:para ratios (more generally, the 
relative populations of the possible nuclear spin states of 
a given molecule or molecular ion) . Section 3 contains the 
results of the calculations of the fractional abundances of 
the ions and of the relative populations of the nuclear 
spin states. Particular attention is paid to the influence of 
the initial value of the ortho:para H 2 ratio. Also, we com- 
pare the results of the calculations, following the grav- 
itational collapse, with the values which are predicted, 
assuming that the medium has time to attain a steady 
state. The possible effects of grain coagulation are consid- 
ered. In Section 4, we summarize our results and make our 
concluding remarks. 

2. The model 

Following the procedures and conclusions of our previous 
paper (Flower et al. 2005), we consider the evolution of a 
contracting sphere of gas and dust and apply the equations 
appropriate to free-fall, homologous, isothermal collapse. 
The timescale for free-fall collapse is 



Tff 



3tt 



32G/90 



(1) 



where po is the initial mass per unit volume. For an ini- 
tial gas density of rin = n(H) + 2n(H 2 ) = 10 4 cm~ 3 , the 
free-fall time is of the order of 10 5 yr. The (constant) 
temperature was T = 10 K. These values of tih and T are 
typical of those prevailing in dark molecular clouds, from 
which protostars form, through gravitational contraction. 

The initial composition was taken to be that of a 
molecular gas in static equilibrium, with the elemental 
abundances and depletions as specified in table 1 of Flower 
et al. (2005). The value of the dust:gas mass ratio was 
0.0094. The chemical reaction set was a composite of our 
previous studies (Walmslcy et al. 2004, Flower et al. 2004, 
2005). It comprises reactions involving species containing 
H, He, C, N, O and S and distinguishes the nuclear spin 
states of H 3 and its deuterated isotopes, in addition to 
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those of and H 2 (see Section |2~T1 below 1 ) . Following the 2.1.1. H2 



discussion in Section 4.4 of Flower et al. (2005), we have 
continued to adopt a grain sticking coefficient 5(N) = 0.1 
for atomic nitrogen, and, with a view to internal consis- 
tency, we have taken 5(0) = 0.1 also. For all other atomic 
and molecular species, S — 1.0. The rate of cosmic ray ion- 
ization of hydrogen was taken to be C = 1 x 10~ 17 s -1 , as 
in our previous paper. 

The initial grain radius was taken to be a g = 0.05 /xm. 
As shown in our previous paper (Flower et al. 2005), this 
value of a g yields approximately the same grain opacity, 
n g 7rag (where n g is the grain number density) , as the grain 
size distribution of Mathis et al. (1977), with limits of 
0.01 < a g < 0.3 fxm. 

The validity of assuming static equilibrium, particu- 
larly for the orthoipara H2 ratio, will be considered be- 
low. In view of the importance of the orthoipara H2 ratio 
- which becomes apparent from the results presented in 
Section [5] - we have varied its initial value between the 
statistical ratio of 3 and the ratio of 3.6 x 10 -7 , which 
corresponds to Boltzmann equilibrium [Equ. J3J below] at 
T=10K. 



2.1. Nuclear spin statistics 

In the calculations reported below, the abundances not 
only of the chemical species but also of their individual 
nuclear spin states have been determined, where applica- 
ble. The cases of interest here involve molecules or molec- 
ular ions comprising two or more protons, such as H2. 
Fcrmi-Dirac statistics require that the nuclear wave func- 
tion should be asymmetric under exchange of identical 
protons, and this restriction associates states of given nu- 
clear spin symmetry with states of appropriate rotational 
symmetry. Thus, in the case of H2, states with total nu- 
clear spin 1 = (para-H2) are associated with rotational 
states with even values of the rotational quantum num- 
ber, J; states with 1=1 (ortho-H2) are associated with 
odd J. Transitions between states of differing nuclear spin, 
1 = and I = 1 in this example, are induced by proton- 
exchange reactions with H + and H3 . [Although HCO + 
is initially the most abundant ion which might undergo 
proton exchange with H2, the measurement by Huntress 
(1977) of the rate coefficient at T = 300 K for the reac- 
tion HCO + (D 2 , HD)DCO + suggests that the analogous 
reaction with H2 is probably too slow to be significant.] 
In view of the importance of H2 in the chemistry of molec- 
ular clouds, it might be anticipated that its orthoipara ra- 
tio would influence the analogous ratios in other molecules 
and molecular ions. 

In order to provide a framework for the interpretation 
of the numerical results which will be presented in the 
following Section, we shall establish first the values of the 
relative populations of the nuclear spin states of a number 
of key species, in static equilibrium, on the basis of an 
analysis which is necessarily statistical in nature. 



We assume that ortho- and para-H2 are formed (on 
grains) in the statistical ratio of 3:1, i.e. the ratio of the 
corresponding values of (21+ 1). Subsequent radiative cas- 
cade leads to the J = 1 (7 = 1) and J = (1 = 0) levels 
of the ground vibrational state, which are separated by 
170.5 K. The rates of cosmic ray dissociation and ioniza- 
tion of H2 are taken to be independent of the total nu- 
clear spin, I. At the low temperature which we consider 
(T = 10 K), proton-exchange reactions convert ortho-H2 
to para-H2; the reverse reaction is insufficiently rapid at 
low T for the orthoipara ratio to attain its Boltzmann 
value, 



n(J = 1) 
n(J = 0) 



9 exp 



-170.5 



(2) 



where the factor of 9 corresponds to the ratio of the sta- 
tistical weights, (21 + 1)(2J +1), of these ortho (J = 1) 
and para (J = 0) levels. At hh = 10 4 cm~ 3 , the steady- 
state value of n(J = l)/n(J = 0) is 2.7 x 10~ 3 , which is 
much larger than the value of this ratio (3.6 x 10~ 7 ) in 
Boltzmann equilibrium at T = 10 K. In other words, at 
low temperatures, ortho-H2 is formed principally by H— 
atom recombination on the surfaces of grains, rather than 
from para-H2 in the gas phase, through proton-exchange 
reactions with H + and Hjj~ . At a higher temperature of 
T = 30 K, the steady-state value of n(J = l)/n(J = 0) 
(3.3 x 10 -2 ) is much closer to its value in Boltzmann equi- 
librium at this temperature (3.1 x 10~ 2 ; cf. Fig. QJ, as 
would be expected. 

2.1.2. H+ 

H^" is formed through cosmic ray ionization of H2, and 
hence the rates of production of ortho- and para-H^ are 
proportional to the number densities of ortho- and para- 
H2, respectively. Both ortho- and para-Hj are destroyed 
by para-H2, which is the much more abundant of the two 
forms of H2, at equal rates, yielding H^". It follows that 
the orthoipara ratio of is essentially the same as that 
of H 2 . 



2.1.3. H+ 

Para-H^" forms principally from para-H^" and para-Bi2, 
which are the more abundant forms of and H2, re- 
spectively, 



H+(p)+H 2 (p) ^H+(p)+H 



(3) 



whereas ortho-H^ is produced from either ortho-H^ and 
para-H2 or para-Hj and ortho-H2, 



H+(o) + H 2 (p)^H+(o) + H 
H+(p)+H 2 (o)^H+(o) + H 



(4) 
(5) 



with the same total rate coefficient (2.1 x 10 9 cm 3 s x ). 
Recalling that the orthoipara ratio of is essentially the 



D.R. Flower et al.: Importance of the orthoipara H2 ratio 



4 

same as that of H 2 , we conclude that the ortho:para H J 
ratio expected from considerations of its formation and 
destruction alone is the same as the ortho:para H 2 ra- 
tio. However, proton-exchange reactions of H;j" with H2 
modify the ortho:para ratio. The proton-exchange re- 
actions of the lowest energy states of ortho- and para-H^" 
with para-H2 are both endoergic, by 137.6 K and 203.4 K, 
respectively, and consequently negligible at T = 10 K. On 
the other hand, the reverse reactions with ortho-H2 

H3 + (o)+H 2 (o)^H 3 + (p)+H 2 (p) (6) 

H 3 ( P ) + + H 2 (o) -> H+ (o) + H 2 (p) (7) 

are exoergic and would yield a ratio of orthoipara H^" of 
4:1, as may be seen from the following argument. 

The statistical weight of the lowest level of para-H^ 
is 2, and that of ortho-H^ is 4. The rate coefficient for a 
reaction is taken proportional to the ratio of the combined 
(by multiplication) statistical weight of the products to 
that of the reactants. Thus, the rate coefficient for the 
reaction (JJJ is taken to be 4 times larger than for reaction 

©■ 

In equilibrium, the population densities of ortho- and 
para-H^" are found to be comparable, i.e. the proton- 
exchange reactions with ortho-H 2 dominate. As a con- 
sequence, the orthoipara ratio depends on the density 
of ortho-H 2 and hence on orthoipara H 2 ratio. 

3. Results 

3.1. Initial abundances 

The elemental abundances and their initial repartition be- 
tween the gas and the solid phases were taken from table 1 
of Flower et al. (2005); the value of the dust:gas mass ratio 
was 0.0094. The chemical composition of the gas, in steady 
state, was determined by running a time-dependent cal- 
culation, at a constant density njj = 10 4 cm -3 and kinetic 
temperature T = 10 K, neglecting any further freeze-out 
on to the grains, until the fractional abundances became 
constant. Under these conditions, the steady-state value 
of the orthoipara H 2 ratio is 2.7 x 10~ 3 . 

Flower & Pineau des Forets (1990) and Pineau des 
Forets et al. (1991) studied the isobaric thermal and 
chemical evolution of interstellar clouds, from initially 
atomic to finally molecular gas which had attained a 
steady state; the values of the gas pressure which they 
considered were similar to the initial value adopted here 
(iihT = 10 5 cm -3 K). Their studies showed that the 
timescale for establishing the steady-state value of the or- 
thoipara H 2 ratio is particularly large, in excess of 10 7 yr. 
At time t « 10 6 yr, the formation of H 2 on grains gives 
rise to a local maximum of the orthoipara ratio, which 
approaches the statistical value of 3:1 associated with 
the formation process. The orthoipara ratio subsequently 
falls towards its steady-state value, which is reached at 
t w 10 7 yr; the steady-state value of the ratio exceeds 
that in Boltzmann equilibrium [Equ. J5J]. Thus, during 



the time interval 10 6 <J t <; 10 7 yr, the orthoipara H 2 
ratio in the molecular gas remains larger than its final, 
steady-state value. 

In Fig. Q] is plotted the ratio of ortho- to para-H 2 as 
a function of time, t, at a constant gas density of nu = 
10 4 cm~ 3 and kinetic temperatures T = 10 K and T = 
30 K; the initial value was taken to be 3, corresponding 
to the statistical ratio associated with the formation of 
H 2 on grains. Fig. ^ shows that the time required for the 
ratio of ortho- to para-H 2 to attain its steady-state value 
(2.7 x 10~ 3 at T = 10 K and 3.3 x 10~ 2 at T = 30 K) is 
approximately 3 x 10 7 yr; this time is comparable with a 
recent determination, by Tassis & Mouschovias (2004), of 
the lifetimes of molecular clouds (« 10 7 yr). [We note that 
Hartmann et al. (2001) estimated the ages of clouds in the 
solar vicinity as being of the order of 10 6 yr.] We conclude 
that it is probable that the orthoipara H 2 ratio does not 
reach its steady-state value before protostellar collapse 
begins. Accordingly, we have varied the initial value of 
the orthoipara H 2 ratio between the statistical value of 3 
and 3.6 x 10" 7 , corresponding to Boltzmann equilibrium 
at T = 10 K. It transpires that the initial value of this 
ratio is critical not only to the populations of the nuclear 
spin states of H^" and its deuterated isotopes but, perhaps 
more significantly, also to the degree of deuteration of 
(cf. Pineau des Forets et al. 1991). 

We see from Fig.^that, for T — 30 K, the orthoipara 
H 2 ratio in steady-state is close to its value in thermo- 
dynamic equilibrium, whereas this is not the case for 
T = 10 K. At low temperatures, the direct formation 
of ortho-H 2 on grains is more rapid than the endoergic 
proton-exchange reactions of H + and with para-H 2 . 
The value of the kinetic temperature, T cr , below which the 
steady-state value of the orthoipara ratio exceeds its value 
in Boltzmann equilibrium may be estimated by equating 
these two rates of formation of ortho-H 2 . Noting that 
ortho- and para-H 2 are assumed to form, on grains, in 
the statistical ratio of 3:1, and equating the total rate of 
formation of H 2 to the rate of its destruction by cosmic 
rays, we obtain T cr w 20 K for the parameters adopted 
in the present calculations. In general, when T <J T cr , the 
initial ortho:para H 2 ratio depends on the thermal history 
of the core. 

3.2. Depletion during collapse 

In Fig. we show the fractional abundances of CO and 
selected isotopes of as functions of the gas density, njj , 
for our 'reference' model of collapse. We note that density 
is related to age through the equation for free-fall collapse. 

The tracer molecule, CO, is already significantly de- 
pleted (through freeze-out) by the time that nn = 3 x 10 4 
cm -3 , whereas the fractional abundances of H 2 D + and D3" 
increase until much higher densities, nn > 3 x 10 5 cm~ 3 , 
corresponding to later times. It is clear that, in this model, 
deuterium fractionation occurs at densities for which CO 
is practically completely depleted (as is N 2 also). In or- 
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Fig. 1. The ratio of abundances of ortho- and para-H2 as 
a function of time, t, from an assumed initial value of 3 
through to the attainment of steady state. A constant gas 
density of riH = 10 4 cm" 3 and temperatures T = 10 K 
and T = 30 K were adopted. The ortho:para ratios in 
Boltzmann equilibrium are shown for both values of T . 



Fig. 2. The fractional abundances of CO and selected iso- 
topes of as functions of the gas density, hh, for a 
constant kinetic temperature, T = 10 K (our 'reference' 
model of collapse). The initial value of the ortho:para H2 
ratio is 2.7 x 10~ 3 . 



der to achieve high levels of deuterium fractionation in 
molecules such as NH3 , it is necessary to minimize the dif- 
ference between the critical densities (or the correspond- 
ing times) at which (i) heavy species, such as CO and N2, 
have frozen out, and (ii) deuterium fractionation, through 
reactions with H2D+ and the other deuterated forms of 
H^", has begun. 

3.3. Steady-state vs. free-fall 

In our previous paper (Flower et al. 2005, fig. 6), we com- 
pared the fractional abundances of H + , H.^ and e~, in 
steady-state and as predicted by the free-fall model, for 
densities 2 x 10 5 < riH < 2 x 10 7 cm~ 3 . The comparison 
was made in the limit of complete heavy element deple- 
tion. At times which are sufficiently large for steady state 
to be attained, essentially all the species containing ele- 
ments heavier than He are frozen on to the grains. The 
agreement between the two calculations was found to be 
good, for the species listed above; but the deuterated iso- 
topes of H3 were not included in the collapse model, and 



neither were the various possible nuclear spin modifica- 
tions. We are now in a position to make such comparisons 
in considerably more detail. However, for reasons which 
will now be given, we have chosen to make the compari- 
son in a somewhat different way from that followed in our 
previous paper. 

The free-fall collapse is initiated at a density jih = 
10 4 cm -3 and for a given repartition of the elements be- 
tween the gas and solid phases (cf. Section 3.1 above). In 
the course of the free-fall collapse, atoms and molecules 
continue to freeze on to the grains, on a timescale which is 
comparable with the free-fall time. Thus, the degree of de- 
pletion of the heavy elements from the gas phase increases 
with time and hence with the density of the medium. 
When comparing the 'collapse' model with the 'steady- 
state' solution, we adopt the composition of the gas at a 
given density on the collapse profile and then compute the 
corresponding steady-state solution, at a time t > 10 8 yr, 
keeping constant both the density and the degree of deple- 
tion of the gas. This procedure ensures that the compari- 
son reflects departures of the free-fall model from steady 
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state for all times (densities), including those which fall 
short of the 'complete depletion' limit. 

3.3.1. Fractional abundances of ions 

First, we compare the abundances of H + and , com- 
puted in steady-state and following the free-fall collapse: 
see Fig. |2K- I n the present calculations, the deuterated 
isotopes of and their nuclear spin modifications have 
been included in the calculations. For the purposes of 
this first comparison, we sum the fractional abundances of 
ortho- and para-H^ . Reactions leading to the deuteration 
of heavy species (those containing elements heavier than 
He, specifically NH3 and N2H+) have not been included 
explicitly in the present model. Nonetheless, qualitative 
conclusions regarding the degrees of deuteration of such 
species can still be drawn, in so far as the deuterated iso- 
topes of H^" are primarily responsible for the deuterium 
fractionation of heavy molecules. We note that the abun- 
dances of the deuterated isotopes of H3 are determined 
mainly by reactions such as H^(HD, H2)H2D + , which 
determine the degree of deuteration of , and not by 
reactions with heavy molecules, such as CO, whose frac- 
tional abundances become rapidly much smaller than that 
of HD, owing to freeze-out. Comparison calculations have 
shown that, at densities rin ;> 10 5 cm' 3 , at which the 
fractional abundances of the deuterated isotopes of 
become significant (cf. Fig.EJ), the reactions with CO have 
no effect on their fractional abundances. 

The results of the two calculations (steady-state and 
free-fall collapse) are in agreement for densities njj J> 
10 6 cm' 3 , for which complete freeze-out has effectively 
occurred. At lower densities, hh ~ 10 5 cm -3 , the steady- 
state solutions underestimate the fractional abundances of 
both H + and H3" - the former by approximately a factor of 
5. Protons are removed by charge transfer reactions with 
neutral species whose ionization potentials are less than 
that of atomic hydrogen. As the heavy neutrals freeze on 
to the grains, the fractional abundance of H + begins to 
increase, and ultimately the steady-state solution merges 
with the results of the collapse model. The timescale for 
establishing ionization equilibrium is small (of the order 
of 1000 yr, for njj = 10 4 cm' 3 ) compared with the equi- 
librium timescale of the chemistry as a whole (Walmsley 
et al. 2004). Consequently, the free electron density cal- 
culated by the collapse model exceeds the steady-state 
value by no more than about 50% over the entire range of 
density, nn- 

The reason that the steady-state solution underesti- 
mates the fractional abundance of Hjj" in the vicinity of 
nn = 10 5 cm' 3 may be seen by referring to Fig. [3Jd, 
where the fractional abundances of the deuterated iso- 
topes of H3 are plotted. In essence, the steady-state so- 
lutions overestimate the degree of deuteration of H^" at 
intermediate densities. The level of agreement between 
the steady-state and free-fall collapse calculations deteri- 
orates as the degree of deuteration of H3 increases, i.e. in 
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Fig. 3. The fractional abundances, relative to nn, of (a) 
H + (full curve), (broken curve) and e _ (chain curve), 
and (b) H2D+ (full curve), D2H+ (long-dashed curve) and 
(short-dashed curve), as functions of nn- The bold 
curves are the predictions of the free-fall model, the light 
curves are the corresponding results in steady-state. 

the sense H2D + , D2H+, D^"; see Fig. [SJd. The differences 
between the steady-state abundances of these species, and 
those calculated following the free-fall collapse, are signif- 
icant in the context of the deuteration of molecules such 
as NH3 , for which deuteron-transfer and proton-deuteron 
exchange reactions with the deuterated isotopes of are 
crucially important. We note that the analogous proton- 
transfer and proton-exchange reactions with H^, H2D + 
and D2H + contribute to ortho-para conversion in NH3. It 
follows that the issues of the degree of deuteration of NH3, 
on the one hand, and the relative population densities of 
its ortho and para nuclear spin states, on the other hand, 
are inter-related (cf. Pineau des Forets et al. 1991). 

3.3.2. Relative population densities of nuclear spin 
states 

In Fig. 0] are compared the relative populations of the two 
lowest nuclear spin states of H+, H 2 D+, D 2 H+ and D+. In 
each case, the ortho:para H2 ratio [specifically, the ratio 
n( J — 1) : n( J — 0)] is plotted, for reference. We see from 



D.R. Flower et al.: Importance of the ortho:para H2 ratio 




10 1 10 s 10* IU' t(f 10 s 10" 10' 



Fig. 4. The relative abundances of the two lowest nuclear 
spin states (in the sense of energetically higher:lower) of 
H 4 ", H 2 D+, D 2 H+ and as functions of 71h, as predicted 
by the free-fall model (bold curves). The corresponding 
values in steady-state are shown as the light curves joining 
the symbols (full circles for H2, full squares for the other 
species). 

Fig. 0] that there are large discrepancies between the pre- 
dictions of the free-fall collapse model and the equivalent 
steady-state calculations. The ortho:para H^" and H 2 D + 
ratios are correlated with that of H 2 , owing to the impor- 
tance of proton exchange with H 2 . Thus, the underestima- 
tion (in steady-state) of the ortho:para H 2 ratio results in 
the ortho:para Hj and H 2 D + ratios being underestimated 
also - by more than an order of magnitude in the case of 
H 2 D + at riH = 10 5 cm -3 . [We note that the ortho:para 
H 2 D + ratio is initially large (at nn = 10 4 cm -3 ), owing 
to proton-exchange reactions with ortho-H 2 , which are 
statistically favourable to the production of ortho-H 2 D + 
(Walmsley et al. 2004).] 

In steady state, the relative value of the population 
densities of the meta and ortho nuclear spin states of Djj" 
is determined principally by deuteron exchange reactions 
with HD, which favour energetically the production of 
ortho-D^" (specifically, the ground state of the ion; the 
first excited state is the lowest state of meta symmetry) . 
However, as noted in Appendix B of Flower et al. (2004), 
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the reaction meta-D^~(ortho-H 2 , HD)ortho-D 2 H + , which 
is endoergic by only 18 K, removes meta-Dj^", even at 
low T. Whilst this reaction is of secondary importance 
in steady state at riu = 10 5 cm -3 , owing to the decrease 
in the ortho:para H 2 ratio with increasing tt-h, it dom- 
inates the removal of meta-D^" in the free-fall model at 
nH = 10 5 cm -3 , for which the ortho:para H 2 ratio is much 
higher than in steady state (see Fig.QJ. As a consequence, 
the meta:ortho D^" ratio is lower at nn — 10 5 cm -3 in the 
free-fall model than in steady state at the same density. 

3.4. Dependence on the initial ortho:para H2 ratio 

The ortho:para H 2 ratio is important in establishing the 
relative populations of the nuclear spin states of a number 
of species, as the discussion in Section [5] should have made 
clear. Furthermore, the initial value of this ratio proves to 
be critical to the degree of deuteration of H^" during free- 
fall collapse. For the purposes of the illustrations which 
follow, grain coagulation has been neglected. 

In Fig. [SI are shown both the fractional abundances 
and the relative populations of the nuclear spin states of 
H^", H 2 D + , D 2 H + and D^". It is apparent from this Figure 
that the fractional abundances of H 2 D + and more highly 
deuterated isotopes of H^" are strongly dependent on the 
initial ortho:para H 2 ratio - whose value is varied between 
3 and 3.6 x 10 -7 (the ratio in Boltzmann equilibrium at 
T = 10 K) - even at densities as high as nn = 10 6 cm" 3 . 
This dependence arises from the reactions 

H 2 D+ (o) + H 2 (o) H+ (p) + HD (8) 

and 

H 2 D+ (o) + H 2 (o) H+ (o) + HD (9) 

which remove directly ortho-H 2 D + and indirectly para- 
H 2 D + , which is converted into ortho-H 2 D + by the 
proton-exchange reaction 

H 2 D+(p)+H 2 (o) ^H 2 D+(o)+H 2 (p) (10) 

which occurs, once again, with ortho-H 2 . Reactions (JSJ 
and © are the reverse of the reactions of ortho- and para- 
H^" with HD, which are responsible for the deuteration 
of H^; as the initial ortho:para H 2 ratio increases, the 
deuteration of Hjj~ is inhibited effectively by these reverse 
reactions. 

From Fig. [H] we see that the populations of the ener- 
getically higher nuclear spin states (i.e. the ortho states) 
of Hjj~ and H 2 D + increase, relative to the lower (para) 
states, as the ortho:para H 2 ratio rises, whereas it is the 
lower states of the mutiply-deuterated isotopes, D 2 H+ 
and D^" , which become relatively more populated. Proton- 
exchange reactions, such as ((TUt and the analogous reac- 
tion involving Hj^, are enabled by the internal excitation 
energy (170 K) of ortho-H 2 , which acts as a reservoir of 
chemical energy. The net effect is the transfer of popu- 
lation to the excited (ortho) states of H 2 D + and H3 . In 
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Fig. 5. The fractional abundances, relative to nn « 
n(H)+2n(H 2 ), of H+, H 2 D+, D 2 H+ and D+ as functions of 
the gas density, nn, at a constant temperature T = 10 K. 
Results are shown for five values of the orthoipara H 2 ra- 
tio, between the statistical value of 3 and 3.6 x 10" 7 , the 
ratio in Boltzmann equilibrium at T = 10 K. The com- 
puted values of the relative populations of nuclear spin 
states of the same species are also shown, in the sense of 
the population of the energetically higher state divided by 
that of the lower. 

the case of H 2 D + , essentially all of the population is in 
the ortho state at low gas densities and high orthoipara 
H 2 ratios. On the other hand, the excited nuclear spin 
states of the multiply-deuterated species, D 2 H + and D^", 
are removed preferentially by ortho-H 2 in the reactions 

D 2 H+(p)+H 2 (o)^H 2 D+(p)+HD (11) 

D+(m) +H 2 (o) -> D 2 H+(o) + HD (12) 

Reaction (|llf) is exothermic, and reaction (|12|) is endother- 
mic, but by only 18 K. (The corresponding reactions of 
ortho-H 2 with ortho-D 2 H + and ortho-D^" are less ener- 
getically favourable, as the reacting ion is in its ground 
state.) 

In Fig. are plotted the fractional abundances of 
ortho-H 2 D + and para-D 2 H + as functions of the gas den- 
sity, nn, in the course of the collapse. Both these species 



Fig. 6. The fractional abundances, relative to tih ~ 
rt(H) + 2n(H 2 ), of ortho-H 2 D + and para-D 2 H + as func- 
tions of the gas density, nn, at a constant temperature 
T = 10 K. Results are shown for five initial (i.e. at 
rtH = 10 4 cm" 3 ) values of the orthoipara H 2 ratio, be- 
tween the statistical value of 3 and 3.6 x 10" 7 , the ratio 
in Boltzmann equilibrium at T = 10 K. The fractional 
abundance of ortho-H 2 D+ in L1544 (Caselli et al. 2003) 
and IRAS 16293E (Vastel et al. 2004), and the fractional 
abundance of para-D 2 H+ in IRAS 16293E (Vastel et al. 
2004), are indicated, as are the densities believed to pre- 
vail in these cores. 



have been observed in protostellar cores, with fractional 
abundances of the order of 10" 10 (Caselli et al. 2003, 
Vastel et al. 2004); the fractional abundances of ortho- 
H 2 D+ in L1544 (Caselli et al. 2003) and IRAS 16293E 
(Vastel et al. 2004), and the fractional abundance of para- 
D 2 H + in IRAS 16293E, deduced from the observations of 
Vastel et al. (2004), are indicated. The orthoipara H 2 ra- 
tio which is inferred from the observations is sensitive to 
the gas density. In the case of IRAS 16293E, if one adopts 
a density of 4 x 10 6 cm" 3 (cf. Lis et al. 2002), all models 
with an initial orthoipara H 2 ratio below 3 are consis- 
tent with the observations. For the lower density of L1544 
(6 x 10 5 < nn < 2 x 10 6 cm -3 ), the data suggest an initial 
orthoipara ratio which is within an order of magnitude of 
its steady-state value of 2.7 x 10" 3 at uh = 10 4 cm" 3 . 
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More definitive conclusions must await not only further 
observations but also the introduction of temperature and 
density gradients into the models of the cores. 



3.5. Influence of coagulation 

The coagulation of grains in the course of the collapse was 
considered in our previous paper (Flower et al. 2005). The 
critical velocity, below which two colliding grains of radius 
a g are assumed to coagulate, was taken to be 

_ 5 

v clit oc a g 6 

following Chokshi et al. (1993). In the extreme cases which 
we considered, the constant of proportionality was taken 
equal to either or 0.4 (with v CI it in cm s _1 and a g in cm), 
corresponding, respectively, to no coagulation or a criti- 
cal velocity for coagulation which is consistent with the 
measurements of Poppe & Blum (1997). The discussion 
so far has been concerned exclusively with the limit of no 

coagulation; now we shall compare with results obtained 

_ 5 

on including coagulation, with w crit = 0.4 a g 6 ; see Fig.0 
In Figs. EH and b are compared the results of calcu- 
lations, with and without coagulation, of the fractional 
abundances of H + , FFj" and its isotopes, and of the free 
electrons. From these Figures, it is apparent that the main 
effect of coagulation is that, at high densities, H + becomes 
the principal ion; this is because the rate of recombina- 
tion of H + with electrons on the surfaces of (negatively 
charged) grains decreases with decreasing grain surface 
area and hence with the degree of coagulation (cf. Flower 
et al. 2005). To the higher density of H+ there corre- 
spond lower values of the orthoipara ratios of H2, Hjj" 
and H2D + , as may be seen from Figs. 0; and d. In ef- 
fect, the excited (ortho) states are more rapidly depop- 
ulated when the density of H + is higher, and so the or- 
thoipara ratios decrease towards their values in Boltzmann 
equilibrium: 9 exp(— 170.5/T) for the orthoipara H 2 ra- 
tio; 2 exp(— 32.9/T) for the orthoipara Hjj" ratio; and 
9 exp(-86.4/T) for the orthoipara H 2 D+ ratio. 

4. Concluding remarks 

We have studied the population densities of molecules 
and molecular ions during the early stages of protostel- 
lar collapse, including H^" and its deuterated isotopes, in 
their various nuclear spin states. The results of a free-fall 
simulation were compared with the corresponding steady- 
state calculation, for densities in the range 10 4 < nn < 
10 7 cm" 3 . Calculations at this level of detail are required 
because current observations of the deuterated isotopes of 
Hjj~ (H2D + , D2H + ) relate to specific nuclear spin states. 
The simulations show that there can be large discrepan- 
cies between the 'free-fall' and 'steady-state' results dur- 
ing the early stages of collapse. 

In general, the computed population densities of the 
different nuclear spin states of molecules and molecular 
ions differ considerably from their values in Boltzmann 




Fig. 7. The fractional abundances of (a) H;j" and (b) its 
isotopes, as functions of nn, as predicted by the free-fall 
model including (bold curves) and excluding (correspond- 
ing light curves) coagulation; see Section 13.51 The frac- 
tional abundances of H + and of the free electrons are also 
plotted in panel (a). The orthoipara ratios of H2, and 
H 2 D + are displayed in panels (c) and (d). 

equilibrium. A striking and important example is H 2 , 
for which the ratio of population densities n(J = 1) : 
n(J = 0), computed in steady state for hh = 10 4 cm" 3 
and T = 10 K, exceeds its value in Boltzmann equilib- 
rium by 4 orders of magnitude. Furthermore, the value of 
the orthoipara H 2 ratio computed following the isothermal 
free-fall collapse exceeds the steady-state value by more 
than an order of magnitude at riu = 10 cm" 3 . It is this 
overabundance which underlies the importance of ortho- 
H2 in inducing population transfer between the different 
nuclear spin states of molecular ions, such H^" and H2D" 1 ", 
at low temperatures. 

In view of the long timescales associated with the pro- 
cesses involved, it seems unlikely that the orthoipara H2 
ratio will reach its steady-state value prior to the com- 
mencement of gravitational collapse. Present observations 
of the protostellar cores L1544 (Caselli et al. 2003) and 
IRAS 16293E (Vastel et al. 2004) are consistent with this 
statement; but the uncertainties in the observations and 
the limitations of the current model are such that defini- 
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tive conclusions cannot be drawn. Our calculations show 
that values of the ortho:para H 2 ratio higher than in 
steady state lead to a reduction in the degree of deuter- 
ation of Hg . As H 2 D+ and D 2 H + are the main agents of 
deutcration of other species, such as N2H + and NH3, a 
decrease in the level of deuteration of H^" implies lower 
levels of deuteration of these other species. High values of 
the ortho:para H2 ratio will occur earlier in the lifetime 
of the progenitor molecular cloud. Thus, one of the pre- 
dictions of the present study is that protostars forming in 
'young' (t <J 10 6 yr) molecular clouds should not display 
high levels of deuteration, owing to the higher prevailing 
values of the ortho:para H 2 ratio. 
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